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Abstract 

We discuss different methods of calculation of the screened Coulomb interaction U in transi- 
tion metals and compare the so-called constraint local-density approximation (LDA) with the GW 
approach. We clarify that they offer complementary methods of treating the screening and, there- 
fore, should serve for different purposes. In the ah initio GW method, the renormalization of 
bare on-site Coulomb interactions between 'id electrons (being of the order of 20-30 eV) occurs 
mainly through the screening by the same 'id electrons, treated in the random phase approxima- 
tion (RPA). The basic difference of the constraint-LDA method from the GW method is that it 
deals with the neutral processes, where the Coulomb interactions are additionally screened by the 
"excited" electron, since it continues to stay in the system. This is the main channel of screening 
by the itinerant {4sp) electrons, which is especially strong in the case of transition metals and 
missing in the GW approach, although the details of this screening may be affected by additional 
approximations, which typically supplement these two methods. The major drawback of the con- 
ventional constraint-LDA method is that it does not allow to treat the energy-dependence of U, 
while the full GW calculations require heavy computations. We propose a promising approxima- 
tion based on the combination of these two methods. First, we take into account the screening of 
Coulomb interactions in the 3(i-electron-line bands located near the Fermi level by the states from 
the subspace being orthogonal to these bands, using the constraint-LDA methods. The obtained 
interactions are further renormalized within the bands near the Fermi level in RPA. This allows 
the energy-dependent screening by electrons near the Fermi level including the same 3d electrons. 

PACS numbers: 71.10.-w, 71.15.-m, 71.20.Be, 79.60.-1 



I. INTRODUCTION 



The description of electronic structure and properties of strongly correlated systems 
presents a great challenge for ab initio electronic structure calculations. The main complexity 
of the problem is related with the fact that such electronic systems typically bear both local- 
ized and itinerant character, where most of conventional methods do not apply. A canonical 
example is the local- [spin]-density approximation (L[S]DA) in the density- functional theory 
(DFT).i 

The DFT, which is a ground-state theory, is based on the minimization of the total 
energy functional E[p] with respect to the electron density p. In the Kohn-Sham (KS) 
scheme, which is typically employed for practical calculations, this procedure is formulated 
as the self-consistent solution of single-particle KS equations 



defined in terms of eigenf unctions (V'i), eigenvalues and the occupation numbers (/j) of 
KS quasiparticles. 

The LSDA provides an explicit expression for Vks[p]- However, it is based on the ho- 
mogeneous electron gas model, and strictly speaking applicable only for itinerant electron 
compounds. 

The recent progress, which gave rise to such directions as LDA+ Hubbard U (Refs. Q0jj) 
and LDA+DMFT (dynamical mean-field theory) (Refs. Oj^, is based on the idea of parti- 
tioning of electronic states. It implies the validity of the following postulates: 

(1) All solutions of KS equations (Q) in LDA can be divided (by introducing proper 
projection-operators) into two subgroups: zG/, for which LSDA works reasonably well, and 
i&L, for which LSDA encounters serious difficulties and needs to be improved (a typical 
example is the 3d states in transition-metal oxides and some transition metals). 

(2) Two orthogonal subspaces, I and L, are "flexible" in the sense that they can be 
defined for a wider class of electron densities, which can be different from the ground- 
state density in LDA. This allows to "improve" LDA by adding a proper correction AE 



{-V^ + VKs[p])^i[p]=eiiJi[p], 



(1) 



which are combined with the equation for the electron density: 




(2) 
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(generally, an tu- dependent self-energy) to the KS equations, which acts solely in the L- 
subspace but may also affect the /-states through the change of p associated with this AE. 
Thus, in the KS equations, the L- and /-states remain decoupled even after including AE: 
(V'iG/[p]|(~V^+yKs[p]+AS)|'?/'jgL[p])=0. For many applications, the L-states are atomic or 
Wannier-type orbitals. In this case, the solution of the problem in the L-space becomes 
equivalent to the solution of a multi-orbital Hubbard-type model, and the formulation of 
the LDA+f/ approach is basically a mapping of the electronic structure in LDA onto this 
Hubbard model. In the following, by referring to the LDA+?7 we will mean not only the 
static version of this method, originally proposed in Ref. Q, but also its recent extensions 
designed to treat dynamics of correlated electrons and employing the same idea of partition- 
ing of the electronic states.^ 

(3) All physical interactions, which contribute to AE, can be formally derived from LDA by 
introducing certain constraining fields {5V^xt} in the subspace of L-states of the KS equa- 
tions (i.e., in a way similar to AE). The purpose of including these {5Kxt} is to simulate 
the change of the electron density, 5p, and then to extract parameters of electronic interac- 
tions from the total energy difference E[p+5 p\—E[p\, by doing a mapping onto the Hubbard 
model. The total energy difference is typically evaluated in LDA,'' and the method itself is 
called the constraint-LDA (CLDA)..S-2iOii 

However, despite a more than decade of rather successful history, the central question 
of LDA+f/ is not completely solved and continues to be the subject of various disputes 
and controversies . This question is how to define the parameter of the effective 
Coulomb interaction U . 

To begin with, the Coulomb U is not uniquely defined quantity, as it strongly depends 
on the property for the description of which we want to correct our LDA scheme. One 
possible strategy is the excited-state properties, associated with the complete removal of an 
electron from (or the addition of the new electron to) the system, i.e. the processes which 
are described by Koopman's theorem in Hartree-Fock calculations and which are corrected 
in the GW method by taking into account the relaxation of the wavefunctions onto the 
created electron hole (or a new electron) . However the goal which is typically pursued 
in LDA+f/ is somewhat different. Namely, one would always like to stay as close as it is 
possible to the description of the ground-state properties. The necessary precondition for 
this, which should be taken into account in the definition of the Coulomb U and all other 
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interactions which may contribute to AS is the conservation of the total number of particles. 
In principle, similar strategy can be applied for the analysis of neutral excitations (e.g., by 
considering the tu-dependence of AS), for which the total number of electrons is conserved.— 
The basic difference between these two processes is that the "excited" electron in the second 
case continues to stay in the system and may additionally screen the Coulomb U. This 
screening may also affect the relaxation effects.— 

The purpose of this paper is to clarify several questions related with the definition of the 
Coulomb interaction U in transition metals. We will discuss both the momentum (q) and 
energy {u) dependence of U, corresponding to the response of the Coulomb potential onto 
the site (R) and time (t) dependent perturbation (5V^xt; and present a comparative analysis of 
the existing methods of calculations of this interaction, like CLDA and GW. We will argue 
that, despite a common believe, the GW method does not take into account the major 
effect of screening of the effective Coulomb interaction U between the 3d electrons by the 
(itinerant) Asp electrons, which may also contribute to the q-dependence of U. This channel 
of screening is included in CLDA, although under an additional approximation separating 
the 3d- and 4sp-states, while in the GW approach, its absence can be compensated by an 
appropriate choice of the pseudo-Wannier orbitals, simulating the basis of L-states. On the 
other hand, CLDA is a static approach, which does not take into account the tu-dependence 
of U.-^ We will consider mainly the ferromagnetic (FM) fee Ni, although similar arguments 
can be applied for other metallic compounds. We start with the basic definition of U for 
the systems with the conserving number of particles, which was originally introduced by 
Herring,— and then discuss the connection of this definition with the parameters which 
comes out from CLDA and GW calculations. 

II. HERRING'S DEFINITION AND CLDA 

According to Herring,— the Coulomb U is nothing but the energy cost for moving a L- 
electron between two atoms, located at R and R', and initially populated by nL-^=nL^i=nL 
electrons: 

f^RR' = E[nL^ + 1, n^R/ - 1] - E[nL^, nj^R/]. (3) 
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In DFT, [/rr' can be expressed in terms of the KS eigenvalues, eLYi=dE /driL^, using Slater's 
transition state arguments:— 

t^RR' = ^Lr[^LR + ^, ^LR' - ^] - ^Lr[^LR " ^LR' + (4) 



The final definition 

S^LR 



(5) 



which is typically used in CLDA calculations, is obtained after replacing the finite differ- 
ence between two KS eigenvalues in Eq. (0} by their derivative. The derivative depends 
on the path in the sublattice of occupation numbers along which it is calculated (e.g., 
^LR+'^LR'=const). This dependence has a clear physical meaning and originates from the 
distance-dependence of intersite Coulomb interactions, which contribute to the screening 
of J/rr'. In the reciprocal (Fourier) space, this distance-dependence gives rise to the q- 
dependence of U . 

Owing to the existence of the second subsystem, /, the reaction Q may compete with 
another one 

U = E[nLR + l,n/R - l,nLR/ - l,n/R/ + 1] - ^[^lr, ^/r, n^R/, n/R/], (6) 

corresponding to independent "charge transfer" excitations at the sites R and R'.— It can 
be also presented in the form (0), but with the different constraint imposed on the numbers 
of L- and /-electrons: nLii+nj-[i=const. Generally, the definitions Q and (jHI) will yield two 
different interaction parameters. Since in the charge-transfer scenario any change of n^R is 
totally screened by the change of n/R located at the same site, the interaction (0) does not 
depend on R. 

In reality, both processes coexist and the proper interaction parameter is given by the 
following equation 

Urr' = E[nLn + 1, n/R - S, ulr' - 1, n/R/ + 5] - E[nLR, n/R, ulr', n/R/], 

where the amount of charge S redistributed between two subsystems is determined varia- 
tionally to minimize J/rr'. In the CLDA scheme, it is convenient to work in the reciprocal 
(Fourier) space and calculate f/q as the response to the q-dependent constraining field 

5Kxt(q,R) = ^LCOsq-R, (7) 



acting in the subspace of L-states under the general condition of conservation of the total 
number of particles. The results of these calculations will strongly depend on how well L- 
electrons are screened by the J-ones. In the case of perfect (100%) screening, the reaction 
(jni will dominate, and the parameter U will not depend on q. If the screening is not 
perfect (e.g., the change of the number of Srf electrons in the transition metals is screened 
to only about 50% by the Asp electrons at the same atom - Ref. llOl ). it is reasonable to 
expect strong q-dependence of the effective [/, because two different channels of screening, 
given by Eqs. Q and (jH)), will work in a different way for different q's. Since the excess 
(or deficiency) of L-electrons caused by a uniform shift of the external potential (5V^xt can 
be only compensated from the system of /-electrons, the "charge transfer" mechanism (jH)) 
will always dominate for small q. The mechanism Q becomes increasingly important near 
the Brillouin zone (BZ) boundary, and will generally compete with the "charge transfer" 
excitations depending on the distribution of the /-electron density.— 



III. THE GW METHOD 
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It was recently suggested by several authors (e.g., in Refs. 
Coulomb U in the LDA-|-f/ approach can be replaced by the screened Coulomb interaction 
W taken from the ah initio GW method. The latter is calculated in the random phase 
approximation (RPA):i^^*iSiii 



W{u) 



uP{uj) 



u. 



(8) 



We adopt the orthogonal atomic-like basis of linear-muffin-tin orbitals (LMTO) 
which specifies all matrix notations in Eq. (jH}. For example, the matrix of bare Coulomb 
interactions e^/|r— r'| has the form {a(3\u\'~^5)=e^ j dr J rfr'x* (r)x^(r')|r— r'|"-'^X^(r)x5(r'), 
and all other matrices are defined in a similar way. The diagonal part of u for the 3d states 
is totally specified by three radial Slater's integrals: F^, F"^, and F^. In the following we will 
identify F^ with the parameter of bare Coulomb interaction, which has the same meaning 
as the Coulomb U after taking into account all screening effects. F'^ and F'^ describe non- 
spherical interactions, responsible for Hund's rule. 

The first advantage of RPA is that it allows to handle the tu- dependence of W, which 
comes from the w-dependence of the polarization matrix P. The most common approx- 



6 



imation for P, which is feasible for ab initio GW calculations, is that of non-interacting 
quasiparticles:— 

PGwr,r,^ => — — — , (9) 

^ u - Sj + ei + i6{fi- fj) 

which is typically evaluated starting with the electronic structure in LSDA (the spin indices 
are already included in the definition of i and j). Generally speaking, the use of -Pgw is an 
additional approximation, which yields a new interaction W^gw- At this stage, it is not clear 
whether it has the same meaning as the effective U derived from CLDA and whether Eq. (jH)) 
includes all necessary channels of screening. It may also include some other effects, which 
should be excluded from the final definition of f/, in order to avoid the double-counting. One 
is the self-screening arising from local (on-site) interactions between the localized electrons. 
These interactions are not accurately treated in RPA.^^ Therefore, the basic idea is to exclude 
these effects from the definition of PVgw aiid to resort this part to the interaction term of the 
Hubbard modeL^^ In this respect, the second important property of RPA is that it allows 
to easily partition different contribution to P and W . If P=Pi+p2 and Wi is the solution 
of Eq. (jHI) for P=Pi, the total W can be obtained from the same equation after substitution 
P^p2 and u-^Wi in Eq. (jH)). For example, if p2=pLL is the part of -Pgw, which includes all 
possible transitions between the localized states, and Pi=Pr is the rest of the polarization, 
the matrix Wr corresponding to Pr, can be used as the interaction part of the Hubbard 
model.i^ 



A. The GW story for fee Ni 

The ferromagnetic fee Ni is the most notorious example where LSDA encounters serious 
difficulties, especially for description of spectroscopic properties. There are three major 
problems:— (i) the bandwidth is too large (overestimated by ~30%); (ii) the exchange 
splitting is too large (overestimated by ~50%); (ii) the absence of the 6 eV satellite. The 
ab initio GW approach corrects only the bandwidth (although with certain tendency to 
overcorrect), whereas the other two problems remain even in GW.-^^^ Therefore, before 
doing any extensions on the basis of GW method, it is very important to have a clear idea 
about its limitations. In this section we would like to clarify several confusing statements 
about screening of W in GW. We argue that the main results of the ab initio GW method 
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can be explained, even quantitatively, by retaining, instead of the full matrix u in Eq. (jH}), 
only the site-diagonal block ull of bare Coulomb interactions between 3d electrons, in the 
atomic-like LMTO basis set. An intuitive reason behind this observation is the form of 
polarization matrix Q, which can interact only with exchange matrix elements. The latter 
are small unless they are calculated between orbitals of the same type, corresponding to 
the self-interaction. The values of radial Slater's integrals calculated in the basis of atomic 
3d orbitals are F°=24.9, F^=ll.l, and F^=6.8 eV, respectively. All other interactions are 
considerably smaller. Hence, it seems to be reasonable to adopt the limit uil^oo, which 
automatically picks up in Eq. (jH} only those matrix elements which are projected onto the 
atomic 3d orbitals, in the LMTO representation. In this sense the ab initio GW method for 
transition metals can be regarded as the RPA solution of the Hubbard model with the bare 
on-site interactions between 3d electrons defined in the basis of LMTO orbitals. In the GW 
method, these interactions are practically not screened by outer electrons. Note, however, 
that the LMTO basis in the transition metals is generally different from the Wannier basis, 
which should be used for the construction the Hubbard Hamiltonian. As it will become clear 
in Sec. I VI If the Wannier representation has several additional features, which may modify 
conclusions of this section to a certain extent. 

Results of these model GW calculations are shown in Fig. ^ In this case, the energy 
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FIG. 1: Characteristic behavior of site-diagonal element of the screened Coulomb interaction 
VF=(xy xy|l4^p{,=o|xy xy) and the matrix element of the self-energy S=(xy|Sq^Q|xy) between xy 
orbitals of the t2g manifold in the F-point of Brillouin zone obtained in the GW approach with the 
bare Coulomb interactions between 3d electrons in the atomic-like LMTO basis set. Inset shows 
amplified near uj=fi. Matrix elements between Cg orbitals show a similar behavior. 

scale is controlled by the bare interaction F^, which predetermines the asymptotic behavior 
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Reiy(oo) (with W denoting the diagonal matrix element of W) and the position of the 
"plasmon peak" of Imiy(ct;) at ~22 eV , which is related with the sharp increase of ReW^(a;) 
at around 25 eV via the Kramers-Kronig transformation. At small u, the behavior of W{uj) 
is well consistent with the strong coupling regime F"— >oo: namely, W{Lj)r^—P^^{uj), which 
is small (~1.8 eV at uj=0) and does not depend on (though it may depend on F"^ and 
F^). All these features are in a good semi- quantitative agreement with results of GW 
calculations .^^^^i^TiiQ 

The self-energy in GW is given by the convolution of W with the one-particle Green 
function G: 

t{uj) = — [ duj'G{u + uj')W{u'). (10) 
2n J 

Therefore, the cj-dependence of S should incorporate the main features of W^u'). Indeed, 
the low-energy part of S (close to the Fermi energy or the chemical potential fi) is mainly 
controlled by Impl^. Since the main poles of ImW and ImG are well separated on the tu-axis 
(the cu-range of ImG is limited by the 3d bandwidth, ~4.5 eV in LSDA for fee Ni, whereas 
the "plasmon peak" of Iml^ is located only at ~22 eV), one has the following relation: 

1 r°° 

dT./duj\^^^ ^- dulmW{u)/uj^. (11) 
^ Jo 

This yields the renormalization factor Z=[l — (9S/(9co'|^^^]~^~0.5, which readily explains the 
reduction of the 3d bandwidth as well as of the intensity of the valence spectrum in ab initio 
GW calculations (Fig.^.-^ 

Away from the Fermi energy (i.e., for energies \uj\ which are much larger than the 3d 
bandwidth), one has another relation ReS(ci;)~— Reiy(ct;), which readily explains the exis- 
tence of the deep minimum of ReS(cc;) near —30 eV as well as large transfer of the spectral 
weight into this region (shown in the inset of Fig. 12)). Therefore, it is not quite right to say 
that the satellite structure is missing in the ab initio GW approach. It may exist, but only 
in the wrong region of u. 

Thus, even besides RPA, the major problem of the GW description for the transition 
metals is the wrong energy scale, which is controlled by the bare on-site Coulomb interaction 
F^ (~20-30 eV) between the 3d electrons. In summarizing this section we would like to 
stress again the following points: 

(1) The major channel of screening of Coulomb interaction in the GW method for the tran- 
sition metals originates from the 3d^3d transitions in the polarization function calculated 
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FIG. 2: The spectral function A{uj)=-^lm.TrG{u;)sgn{u;-fi) for fee Ni in LSDA and two GW 
schemes with bare electronic interactions and parameters extracted from constraint-LDA. The inset 
shows the satellite structure in A(io) at the F-point of Brillouin zone in the bare-GW approach. 

in the atomic-like LMTO basis set. The screening by the 4sp-electrons is practically absent; 

(2) At small u, the deficiency of the 3d-4sp screening is masked by the strong-coupling 
regime realized in RPA equations for screened Coulomb interaction, which explains a small 
value of W{0) obtained in the GW calculations; 

(3) The main u;-dependence of E and W in GW also comes from the 3d^3d transitions. 
Different conclusions obtained in Refs. are related with the use of different par- 
titioning into what is called the "Sc?" and "non-3(i" (pseudo-) Wannier orbitals.— In the 
light of analysis presented in this section, the strong cu-dependent screening by the "non-Sc?" 



Wannier states obtained in Refs. 



171 means that in reality these states had a substantial 



weight of "3(i" character of the LMTO basis, which mainly contributed to the screening. 
We will return to this problem in Sec. IVIIl 

The next important interaction, which contribute to the screening of F'^ in GW is due 
to transitions between states with the same angular momentum: i.e., 3d^nd {n= 4, 5, 



...) (see also comments in Sec. IV A|) . In the lowest order (non-self-consistent RPA), these 
contributions can be estimated as 



Aiy(a;) ^ {3d3d\u\3dAd)l^PG^{uj,3d Ad) + (higher n), 



(12) 
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where (3(i3(i|'u|3(i4(i)av— 6.1 eV is the spherical part of the exchange integral {3d3d\u\3d4:d) , 
corresponding to Results of these calculations are shown in Fig. El The region of 3d^Ad 

8 

^ 4 
> 

-8 

20 40 60 80 100 

co(eV) 

FIG. 3: The w-dependence of on-site Coulomb interaction associated with the relaxation of the 
3d wavefunctions in the region of the 3(i— >4(i transitions. The 3(i— >5(i transitions have been also 
taken into account. They contribute to the region above 100 eV, which is not shown here. 

transitions strongly overlaps with the "plasmon peak" of ImVF(ti;) (Fig. Q). Therefore, in 
the GW calculations, these two effects are strongly mixed.— li^ji^^ The w-dependence of 
AW will also contribute to the renormalization of the low-energy part of spectrum. In GW, 
this contribution can be estimated using Eq. (jllj) . which yields 9S/9ti;|^^^~0.06. This 
contribution is small and can be neglected. 

IV. GW VERSUS CLDA 

What is missing in the ab initio GW method, and what is the relation between GW and 
CLDA? Let us consider for simplicity the static case, where 5l4xt does not depend on time 
(the generalization to the time-dependent case is rather straightforward). 

Eventually, both methods are designed to treat the response 6p{r) of the charge density 
to the change of the external potential SVext, which can be calculated in the first order of 
the regular perturbation theory. Then, 6Vext will affect both eigenvalues and eigenfunctions 
of the KS equations (P). The corresponding corrections are given by the matrix elements 
{ipilSVextl'ipj) with i=j and ij^j, respectively. If two (or more) eigenvalues are located near 
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the Fermi level, their shift can lead to the repopulation effects when some levels become 
occupied at the expense of the other ones. This is a direct consequence of the conservation 
of the total number of particles, which affects the occupation numbers. Therefore, very 
generally, the total response 6p{r) in metals will consist of two parts, 6p{r)=6ip{r)+62p{r), 
describing the change of the occupation numbers, (5ip(r)=^. 5/j|-?/'j(r)p, and the relaxation 
of the wavefunction, S2p{r)=J2i /t'^|V'i(r)P, respectively. Then, the polarization function P, 
defines as 



will also consist of two parts. Pi and P2, which yield 6ip and 62P after acting on 6Vext- 
Then, it is easy to verify by considering the perturbation-theory expansion for {ipi} with 
respect to 6Vext that the GW approximation corresponds to the choice Pi=0 and P2=Pgw- 
It yields 62p{r), which further induces the new change of the Coulomb (Hartree) potential 
52^(r)=e^/ dr'62p{r')/\r—r'\. By solving this problem self-consistently and taking the 
functional derivative with respect to 62P one obtains the GW expression (jH)) for the screened 
Coulomb interaction VFgw(O). Therefore, it is clear that the ab initio GW method takes into 
account only one part of the total response 6p, describing the relaxation of the wavefunction 
with the fixed occupation numbers. Another contribution, corresponding to the change of 
the occupation numbers (or the charge redistribution near the Fermi level) is totally missing. 

This result can be paraphrased in a different way, which clearly illustrates its connection 
with the definition of orthogonal subspaces, L and J, discussed in the introduction, and the 
partitioning of the polarization function P (Sec. IIII|) . which is used in the definition of the 
Hubbard model . ^^1^^ First, recall that according to the main idea of the LDA+f/ method (see 
postulates 1-3 of the Introduction part), 6Vcxt should be a projector-type operator acting 
in the subspace of the L states. Then, the result of the action of the polarization function 
Pgw=P2, given by Eq. Q, onto this SVext will belong to the same L space. Therefore, the 
projection (5V^xt will generate only that part of the polarization function, which is associated 
with the transitions between localized states {Pll in Sec. IIII|) . Meanwhile, this polarization 
effect should be excluded from the final definition of the parameter U in the Hubbard model 
to avoid the double counting , ^^i^'' However, if Pn is excluded, there will be nothing left 
in the polarization function that can interact with 6Vext and screen the change of the 
electron density in the L-subspace. Therefore, the GW scheme should correspond to the 
bare Coulomb interaction, that is totally consistent with the analysis presented in Sec. IIII Al 




(13) 
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A. Basic Difficulties for Transition Metals 



There is certain ambiguity in the construction of the Hubbard model for the transition 
metals, which is related with the fact that their LDA electronic structure cannot be described 
in terms of fully separated L- and /-states without additional approximations. In this section 
we briefly review two such approximations, which will explain the difference of our point of 
view on the screening of Coulomb interactions in the transition metals from the one proposed 
in Refs.yandbj 

The GW approach employed in Refs. QandH implies that all electronic structure near 
the Fermi level can be described in terms of only five pseudo-Wannier orbitals of predom- 
inantly 3(i-character, which serve as the L-states in the considered model. Generally, such 
L-states are not the same as the LMTO basis functions and take into account the effects of 
hybridization between Be? and 4sp states. An example of such electronic structure, obtained 
after elimination of the 4sp-states near the Fermi level through the downfolding procedure,— 
is shown in Fig. 01 Other possibilities of defining these pseudo-Wannier functions, which 



28l . Then, the remain- 



have been actually used in Refs. QandH are summarized in Ref. 
ing electronic states, which are orthogonal to these pseudo-Wannier orbitals, represent the 
/-states. By the construction, the /-states are expected to be far from the Fermi level. This 
may justify the use of the GW approximation for the screening of Coulomb interactions in 
the 3(i-electron-like bands, formed by the pseudo-Wannier orbitals near the Fermi level, by 
the remote /-states. The parameters of Coulomb interactions, constructed in such a way, 
correspond to the original Herring definition (jH)) in the basis of pseudo-Wannier orbitals. 
Formally, it should also include the charge redistribution effects near the Fermi level. How- 
ever, in this case the charge redistribution goes between pseudo-Wannier orbitals of the same 
(L) type, which constitutes the basis of the Hubbard model. Therefore, the effects of the 
charge redistribution can be taken into account by including the intra- as well as inter-site 
Coulomb interactions in the Hubbard Hamiltonian. The latter can be evaluated in the GW 
approach, provided that the relaxation effect are not very sensitive to whether the excited 
electron is placed on another L-orbital of the same system, or completely removed from it, 
like in the GW method. 

The model employed in CLDA calculations is obtained after neglecting the hybridization 
between 3(i- and 4sp-states (the so-called canonical-bands approximation-Ref . |25|) . It con- 
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FIG. 4: Two approximate views on the electronic structure of (paramagnetic) fee Ni underlying 
different schemes of calculation of the screened Coulomb interaction. The original LDA bands are 
shown by light color. The GW calculations are based on the model 'a', which implies that all 
electronic structure near the Fermi level (located at zero) can be described in terms of only five 
pseudo-Wannier orbitals of predominantly 3d-character, simulating the L-states. The dark bands 
show an example of such electronic structure obtained after elimination of 4sp-states through the 
downfolding procedure)^ The remaining electronic states are the /-states, which are orthogonal 
to the pseudo-Wannier orbitals and allowed to screen the Coulomb interactions in these bands. 
The screening is treated in RPA. The model 'b', which is used in constraint-LDA calculations, is 
obtained after neglect ing the hybridization between 3d- and 4sp-states (the so-called canonical- 
bands approximation 25). It consists of the 3d band (representing the L-states and shown by 
dark color), embedded into the free-electron- like 4sp-band (representing the /-states and shown by 
dash-dotted line). The coexistence of two different groups of states near the Fermi level gives rise 
to the charge redistribution, which contribute to the screening of Coulomb interactions in the 3d 
band. 

sists of the pure Sd-hand, located near the Fermi level and representing the L-states of the 
model, which is embedded into the free-electron-like 4sj9-band, representing the /-states.— 
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Formally, these bands are decoupled and the free-electron-line 4sp-band can be eliminated 
from the basis in the process of construction of the Hubbard Hamiltonian. However, in this 
case the definition of the screened Coulomb interaction in the 3d band should take into ac- 
count the processes corresponding to redistribution of electrons between 3d- and Asp-hand at 
the low-energy cost, which is traced back to Herring's scenario of screening in the transition 
metals,— and which is missing in the GW method. 

However, we would like to emphasize again that both considered models are approxima- 
tions to the real electronic structure of fee Ni. Even in the first case (model 'a' in Fig. |^, 
the free-electron-like 4sp-band lies near the Fermi level (especially around L-point of the 
Brillouin zone). Therefore, the charge redistribution effects are expected to play some role 
even in the basis of Wannier orbitals. On the other hand, because of strong hybridiza- 
tion between 3d- and 4sp-states in the transition metals, there is a substantial difference 
of electronic structure used in CLDA calculations (model 'b' in Fig. from the real LDA 
electronic structure of fee Ni. Strictly speaking, all partial contributions to the screening of 
Coulomb interactions, which we will consider in the next section, will be evaluated for this 
particular model of the electronic structure. The values of these parameters can be revised 
to a certain extent after taking into account the hybridization between 3d- and 4sp-states. 
For example, with the better choice of the Wannier basis for the five 3(i-electron-line bands 
in the model 'b' one could possibly incorporate the main effects of the model 'a' and merge 
these two approaches. 



V. CLDA FOR TRANSITION METALS 



How important are the relaxation of the wavefunctions and the change of the occupation 
numbers in the definition of the Coulomb interaction ?7? For the transition metals, both 
contributions can be easily evaluated in CLDA. For these purposes it is convenient to use 
the Hellman-Feinman theorem, which relates the static U with the expectation value of the 
KS potential Vks='^h+Vxc:— 

U = {3d\^\3d). 

Then, the exchange-correlation (XC) part is small. (5Vh can be expressed through Sp. Hence, 
the CLDA scheme provides the self-consistent solution for 6p associated with the change of 
the number of 3d electrons, dn^d- The latter is controlled by 6Vext- Therefore, the procedure 



15 



is totally equivalent to the calculation of the polarization function P and the screened 
Coulomb interaction for uj=0. 

A. Conventions 

We use rather standard set up for the CLDA calculations. Namely, the 3d band of Ni 
should be well separated from the rest of the spectrum (otherwise, the LDA+f/ strategy 
discussed in the Introduction does not apply). For fee Ni this is not the case. However, 
this property can be enforced by using the canonical bands approximation in the LMTO 
method.— We employ even cruder approximation and replace the 3d band by the atomic 
3d levels embedded into the Asp band (in the other words, we switch off the hybridization 
between 3d orbitals located at different atomic sites as well as the 3d and Asp states).— 
Then, each 3d orbital can be assigned to a single atomic site. By changing the number 
of 3(i-electrons at different atomic sites {R} in supercell calculations, one can mimic the 
q-dependence of the external potential ((Zj). Other atomic population (of the Asp states) are 
allowed to relax self-consistently onto each change of the number of 3d electrons. Hence, 
the contribution of the charge-transfer excitation © to the screening of U is unambiguously 
defined by the form of the external potential and details of the electronic structure of the 
Asp states. Some aspects of treating the 3d states beyond the atomic approximation will be 
considered in Sec. IVllI 

The LMTO method is supplemented with an additional spherical approximation for 
VKs(r) inside atomic spheres, which bars small exchange interactions between 3d and Asp 
electrons from the screening of U. By paraphrasing this statement in terms of the polar- 
ization function in the GW method, the spherical approximation for VKs(r) in the CLDA 
calculations is equivalent to retaining in Pqw only those contributions which are associated 
with transitions between states with the same angular momentum (e.g., 3d—^Ad, etc.). 

B. Screened Coulomb Interaction in the F-point 

First, we evaluate the pure effect associated with the change of the occupation numbers, 
without relaxation of the wavefunctions. This mechanism is directly related with the con- 
servation of the total number of particles, and simply means that the excess (or deficiency) 
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of the 3d electrons for q=0 is always compensated by the Asp electrons, which participate 
in the screening of 3d interactions. The corresponding contribution to the screening of 
is given by:— 

In transition metals, A^^^U is very large and takes into account more than 70% of screening 
of the bare Coulomb interaction (Table This contribution is missing in the GW 
method. The second largest effect (~25% of the total screening) is caused by relaxation 
of the 3d orbitals onto the change of the Hartree potential associated with the change of 
these occupation numbers {A^'^^U in Table E)). The remaining part of the screening (~5%) 
comes from the relaxation of other orbitals (including the core ones) and the change of the 
XC potential. In principle, the relaxation effects should be taken into account by the GW 
calculations. However, this procedure strongly depends on the way how it is implemented. 
For example, the CLDA approach is based on a direct solution of KS equations supplemented 

TABLE I: Partial contributions to the screening of the 3d interactions in the F-point extracted 
from constraint-LDA calculations (in eV): (1) bare Coulomb integral F^, (2) the screening of F^ 
by the Asp electrons associated with the change of occupation numbers, without relaxation of the 
wavefunctions (A^^^F^), (3) the additional screening of F'^ associated with relaxation of the 3d 
orbitals (A^^^F*^), and (4) the total value of U obtained in CLDA calculations. 



compound 


po 


A«F0 


A(2) F° 


U 


bcc Fe 


22.2 


-13.6 


-3.5 


4.5 


fee Ni 


24.9 


-14.2 


-5.2 


5.0 



with a flexible atomic basis set, like in the LMTO method.-^ Then, the change of F^ caused 
by relaxation of the 3d orbitals can be easily evaluated as^ 

2 dn^d 

Since n^d is large in the fee Ni, this contribution is also large. The situation can be different 
in the GW scheme, based on the perturbation theory expansion, which requires a large basis 
set.— For example, in order to describe properly the same relaxation of the 3d wavefunctions, 
the polarization Pgw should explicitly include the excitation from the occupied 3d to the 
unoccupied Ad (and probably higher) states.— 
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C. q-dependence of Coulomb U 



Since the change of the number of 3d electrons in transition metals is not totally screened 
by the Asp electrons at the same atomic site it is reasonable to expect an appreciable q- 
dependence of the effective U. Results of CLDA calculations for the high-symmetry points 
of the Brillouin zone are summarized in Table |n] The effective U appears to be small in the 
F-point due to the perfect screening by the Asp electrons. At the Brillouin zone boundary 
this channel of screening is strongly suppressed that is reflected in the larger values of the 
Coulomb U. The screening by intersite Coulomb interactions, which takes place in the 
X-point of the BZ, is substantially weaker and cannot fully compensate the lack of the 4sp- 
screening. In the L-point of the BZ for the fee lattice, the modulation of the 3(i-electron 
density in the CLDA calculations is such that the number of nearest neighbors with excessive 
and deficient number of 3d electrons is the same. Therefore, the contributions of intersite 
Coulomb interactions to the screening are cancelled out, resulting in the largest value of the 
effective U in this point of the BZ. 

TABLE IL Coulomb interaction U (in eV) for fee Ni in three different points of the Brillouin zone: 



r=(0,0,0), X=(27r,0,0), and L=(7r,7r,7r) 


(in units of 1/a, 


where a is the cubic lattice parameter). 


r 


X 


L 


5.0 


6.8 


7.3 



VI. GW STARTING WITH CLDA 

In this section we discuss some relevance of parameters of effective Coulomb interactions 
extracted from CLDA for the analysis of electronic structure and properties of fee Ni. We 
consider the "renormalized GW approach" , in which, instead of bare Coulomb interactions, 
we use parameters extracted from CLDA. The main difference is that the latter incorporates 
the screening by the 4sp-electrons, including the effects of charge redistribution beyond the 
GW approximation. This strategy can be well justified within RPA, because it allows to 
partition the polarization function and treat the screening effects in two steps: 
(1) We take into account the screening by "non-Sc?" electrons using CLDA. This yields 
the new ("renormalized") matrix of screened Coulomb interactions ull between the 3d 
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electrons.— As it was discussed in Sec. IV CI the obtained interaction un is q-dependent, 
and this dependence is fully taken into account in our calculations. 

(2) We evaluate the screening caused by 3c?— >3c? transitions in the polarization function Q 
using Eq. (jSJ in which the matrix of bare Coulomb interactions un is replaced by un. This 
yields the new interaction W{uj), which is used in subsequent calculations of the self-energy 
fllO|) . It is reasonable to expect that the main tu-dependence of W will come from the 3(i— s>3(i 
transitions (see closing arguments in Sec. IIII which are taken into account in the second 
step. The screening by "non-3(i" states can be treated as static. 

Results of these calculations are shown in Fig. ^ The main effect of the 4sp-screening, 

12 
10 
8 
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FIG. 5: The same as Fig. ^ but with the parameters of Coulomb interactions extracted from 
CLDA. 

beyond the standard GW approach, is the change of the energy scale, which is now controlled 
by the q-dependent Coulomb interaction f/, being of the order of 5.0-7.3 eV. It change the 
asymptotic behavior ReH^(oo) as well as the position and the intensity of the "plasmon 
peak" of Iml^(ci;), which is shifted to the lower-energies region and becomes substantially 
broader in comparison with the case of bare Coulomb interactions considered in Sec. IIII Al 
On the other hand, the static limit Reiy~1.9 eV is practically not affected by details of 
the 4sp-screening, due to the strong-coupling regime realized in the low-u; region. The HeW 
exhibits a strong cu-dependence at around 7 eV, which is related with the position of the 
plasmon peak of \mW{uj). All these features are well reflected in the behavior of 'L{uj). 

The main effect of the 4sp-screening onto the spectral function in RPA consists in some- 
what milder reduction of the bandwidth, which is also related with the spectral weight 
transfer (Fig. |2I): the new renormalization factor is Z~0.7 against Z~0.5 obtained with 
bare Coulomb interactions. However, the exchange splitting does not change and the 6 eV 
satellite structure does not emerge. 
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VII. SUMMARY AND REMAINING QUESTIONS 



We have considered several mechanisms of screening of the bare Coulomb interactions 
between 3d electrons in transition metals. We have also discussed different methods of 
calculations of the screened Coulomb interactions. Our main results can be summarized as 
follows. 

(1) The processes which mainly contribute to the screening of Coulomb interactions between 
3d electrons are essentially local, meaning that the on-site Coulomb interactions are most 
efficiently screened by the 3d and Asp electrons located at the same site .^i-^^i-^^ The most 
efficient mechanism of screening is basically the self-screening by the same 3d electrons, 
evaluated in some appropriate atomic-like basis set, like that of the LMTO method employed 
in the present work. The cu-dependence of the effective Coulomb interaction U also originates 
mainly from the self-screening. 

(2) We have clarified a fundamental difference between constraint-LDA and GW methods 
in calculating the effective Coulomb interaction U. The GW approximation does not take 
into account a screening of the on-site Coulomb interactions by the itinerant Asp electrons, 
taking place via redistribution of electrons between 3d and Asp bands. 

In a number of cases, the GW approach may be justified by using Wannier basis functions, 
representing the bands near the Fermi level. If these bands are well isolated from the other 
bands, the redistribution of electrons between Wannier orbitals for the bands near the Fermi 
level and those far from the Fermi level must be negligible. Then, the remote bands can 
participate in the screening of Coulomb interactions in the "near-Fermi-level bands" only 
via virtual excitations, which can be treated on the RPA level. 

However, in the case of Ni, such separation of bands is not complete, and it is essential to 
consider additional mechanisms of screening beyond the GW approximation. In the present 
work, the 4sp-screening is automatically taken into account in the CLDA approach, which 
is complementary to the GW method. Due to the strong-coupling regime realized in RPA 
equations for the screened Coulomb interaction, the static limit appears to be insensitive to 
the details of the 4sp-screening. However, from the viewpoint of the present approach, the 
4sp-screening becomes increasingly important at finite u and controls both the asymptotic 
behavior and the position of the plasmon peak of the screened Coulomb interaction in RPA. 
The latter effect can be especially important as it predetermines the position of the satellite 
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structure. 

Finally, we would like to make several comments about implication of the parameters of 
screened Coulomb interaction obtained in our work for the description of electronic structure 
and properties of transition metals. We will also discuss some future directions and make a 
comparison with already existing works. 

(1) Our results clearly show that RPA is not an adequate approximation for the electronic 
structure of fee Ni. Even after taking into account the additional screening of the 3d-3d 
interactions by the itinerant Asp electrons, beyond the GW approximation, and the q- 
dependence of the effective U, we obtain only a partial agreement with the experimental 
data. Namely, only the bandwidth is corrected in this "renormalized GW approach", in 
a better agreement with the experimental data. However, there is only a tiny change of 
the spectral weight around 6 eV (Fig. |21), i.e. in the region where the satellite structure 
is expected experimentally. Even assuming that our parameters of Coulomb interactions 
may be still overestimated (due to the reasons which will be discussed below), and the 
satellite peak can emerge for some smaller values of f/,^- one can hardly expect the strong 
spin-dependence of this satellite structure as well as the reduction of the exchange splitting, 
which are clearly seen in the experiment,— on the level of RPA calculations. Therefore, it 
is essential to go beyond. 

(2) Even beyond LDA, do the parameters of screened Coulomb interaction [/~5. 0-7.3 eV, 
obtained in the atomic approximation, provide a coherent description for the electronic 
structure and properties of fee Ni? Probably, this is still an open question because so far 
not all of the possibilities in this direction have been fully investigated. One new aspect 
suggested by our calculations is the q-dependence of the effective U. On the other hand, 
all previous calculations suggest that the Coulomb interaction of the order of 5.0-7.3 eV is 
probably too large. For example, the value of U, which provide a coherent description for a 
number of electronic and magnetic properties of fee Ni on the level of DMFT calculations 
is about 3 eV,-^ which is well consistent with the previous estimates based on the t-matrix 
approach.i^ Therefore, it is reasonable to ask if there is an additional mechanism, which 
further reduces the effective U from 5.0-7.3 eV till 3.0 eV? One possibility lies in the atomic 
approximation which neglects the hybridization effects between 3d and Asp states, and which 
is rather crude approximation for the transition metals.— The hybridization will generally 
mix the states of the 3d and Asp character, and therefore will affect the form of the Wannier 
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orbitals constructed from the atomic wavef unctions. Since the 3d, 4s, and Ap states belong to 
different representations of point group of the cubic symmetry, they cannot mix at the same 
site. However, the 4s (or 4p) orbital can have tails of the 3d character at the neighboring 
sites (and vice versa). These tails will additionally screen the Coulomb interactions between 
the (nominally) 3d electrons. The screening is expected to be very efficient because it 
operates between orbitals of the same {3d) type. It should explain further reduction of the 
static U obtained in the atomic approximation. Another feature of this screening is the lj- 
dependence of the effective U, which comes from the 3d^3d transitions in the polarization 
function (namely between tails of the 4sp-orbitals and the heads of the wavefunctions of 
the 3d character). In RPA, this cu-dependence is directly related with the static limit of 
screening via the Kramers- Kronig transformation.— We believe that the screening by the 
tails of the Wannier functions was the main phy sical mechanism underlying the calculations 
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although this idea has not been clearly spelled out before. The effect of charge redistribution 
between different states located near the Fermi level, which is not taken into account in the 
GW approximation, is also expected to be smaller with the proper choice of the Wannier 
orbitals. 

Another problem is that the 3d and Asp bands are strongly mixed in the case of pure 
transition metals. Therefore, the construction of the separate Wannier functions of the "3(i" 
and "non-3(i" type will always suffer from some ambiguities.— In this sense, the transition- 
metal oxides, whose physical properties are mainly predetermined by the behavior of a 
limited number of 3d bands, located near the Fermi level and well separated from the rest of 
the spectrum, are much more interesting systems for the exploration of the idea of screening 
of Coulomb interactions, formulated on the Wannier basis. For example, based on the above 
argument, one can expect a very efficient screening of Coulomb interactions in the 3d band 
by the Wannier states constructed from the oxygen 2p orbitals, which have appreciable tails 
of the 3d character at the transition-metal sites. The first attempt to consider this screening 
have been undertaken in Ref. on the basis of constraint-LDA method. Similar scheme 
can be formulated within RPA, which takes into account the ^-dependence of the screened 
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Coulomb interaction U. This work is currently in progress. 
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